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Abstract. I review recent developments in the field of relativistic hydrodynamics and 
its application to the bulk dynamics in heavy-ion collisions at the Relativistic Heavy- 
Ion Collider (RHIC) and the Large Hadron Collider (LHC). In particular, I report on 
progress in going beyond second order relativistic viscous hydrodynamics for conformal 
fluids, including temperature dependent shear viscosity to entropy density ratios, 
as well as coupling hydrodynamic calculations to microscopic hadronic rescattering 
models. I describe event-by-event hydrodynamic simulations and their ability to 
compute higher harmonic flow coefficients. Combined comparisons of all harmonics to 
recent experimental data from both RHIC and LHC will potentially allow to determine 
the desired details of the initial state and the medium properties of the quark-gluon 
plasma produced in heavy-ion collisions. 



1. Introduction 

The observed large elliptic flow at RHIC and recently at the LHC is one of the most 
striking observations in heavy-ion collision experiments. This asymmetry of particle 
production in the transverse plane of the collision is interpreted as the hydrodynamic 
response to the initial geometry. The applicability of hydrodynamics demands a short 
mean free path with respect to the system size. Therefore it is concluded that the 
created quark-gluon plasma is strongly coupled and behaves like a nearly perfect fluid. 
Recently interest has cascaded to all higher harmonics, including odd ones such as the 
triangular flow coefficient fs, which are non-zero in single events. A lot of progress has 
been made in extracting medium properties from hydrodynamic calculations, with the 
largest uncertainty emerging from the limited knowledge of the initial conditions. 

2. Theoretical framework 

The current standard for the viscous hydrodynamic description of relativistic heavy-ion 
collisions has been established in [1], where all terms up to second order in gradients 
for a conformal fluid have been derived. Additional terms for non-conformal fluids with 
non-zero bulk viscosity have been derived in [2]. 

In the ideal case, the evolution of the system created in relativistic heavy-ion 
collisions is described by the following 5 conservation equations 

5.^r = o, « = o, (1) 
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where Tl^^ is the energy-momentum tensor and J'^ is the net baryon current. These are 
usually re-expressed using the time-like flow four-vector as 

Tt/ = {e + V)u^u''-Vg^\ J's = Pb^^ , (2) 

where e is the energy density, V is the pressure, ps is the baryon density and 
gt^'^ = diag(l, —1, —1, —1) is the metric tensor. The equations are then closed by adding 
the equilibrium equation of state 

V = V{e,pB) (3) 
local constraint on the variables. 

In the flrst-order, or Navier-Stokes formalism for viscous hydrodynamics, the stress- 
energy tensor is decomposed into Tl'^^ = Tj^^ + S^^ , where T-^"^ is given in Eq. ([2]) and 
the viscous part of the stress energy tensor is given by 

S''" = J] [s/^u" + S/'u^ - ^A'^'^V^n"^ , (4) 

where A'^'^ = g^" — u'^u^ is the local 3-metric and = A^'^di, is the local spatial 
derivative. Note that S'^" is transverse with respect to the flow velocity since A^^'n^ = 
and u^Uiy = 1. Hence, u'^ is also an eigenvector of the whole stress-energy tensor with 
the same eigenvalue e. r] is the shear viscosity of the medium. 

The Navier-Stokes form is conceptually simple but introduces unphysical super- 
luminal signals that lead to numerical instabilities. 

The second-order Israel- Steward formalism [31 HI [5] avoids this super-luminal 
propagation, as does the more recent approach in [6]. In the Israel-Stewart formalism 
for a conformal fluid, derived in [1], the stress-energy tensor is decomposed as T'^" = 
'^id + '^^'^ ■ "^^^ evolution equations are 

9^r^" = o, (5) 

A^A^'^S^tt"^ = - - (tt^'^ - S^") - ^TT^^'id^u'^) , (6) 

where we neglected vorticity and terms that turn out to be numerically irrelevant. For 
the role of vorticity in heavy- ion collisions when including fluctuations see [7j. 

Simulations of bulk dynamics in heavy-ion collisions using this formalism have been 
performed in 2+1 dimensions in [HI [9l [101 [HI [121 [13] , and within the equivalent Ottinger- 
Grmela [13] formalism in [15]. Recently, 3-1-1 dimensional viscous calculations have also 
become available [Ml [17] . 

Having established the standard theoretical framework, in the following I will 
discuss recent developments in the fleld of relativistic viscous hydrodynamics and its 
application to heavy-ion collisions. 

3. Equation of state 

Typically, the equation of state ([3]) used in hydrodynamic simulations of heavy-ion 
collisions is determined from lattice QCD calculations combined with a hadron gas 



Flow in heavy-ion collisions - Theory Perspective 



3 



model. A recent parametrization and a detailed comparison of different lattice equations 
of state has been performed in [18]. The different equations of state [191 [201 1211 [221 [T3]. 
which have different results for the trace anomaly and the speed of sound, lead to 
different evolution of the momentum anisotropy when used in an ideal hydrodynamic 
evolution. However, the difference in the final spectra and elliptic flow results turns out 
to be negligible. This leads to the conclusion that one cannot determine details about 
the equation of state from the comparison of hydrodynamic calculations to experimental 
data. On the other hand, using one of the latest lattice equations of state will do 
a good enough job when extracting medium properties such as transport coefficients 
from hadronic observables. Thermal photon production, however, is potentially more 
sensitive to the equation of state, because photons are produced throughout the 
evolution and hence probe the dynamics of the system more directly [23]. Another 
relevant aspect is the inclusion of chemical freeze-out in the equation of state (and the 
freeze-out procedure) to reproduce the correct final particle ratios [211 [2S1 126] . 

4. Bulk viscosity and limits of second order viscous hydrodynamics 

QCD is a non-conformal theory and contains a finite bulk viscosity. In [27] the effect 
of bulk viscosity on elliptic flow and in [281 129] the combined effect of shear and bulk 
viscosity has been studied. Bulk viscosity is expected to peak (possibly along with 
its relaxation time) around the critical temperature Tc, where the system may develop 
large correlation lengths [301 Ell ESI ES]. In [29], within the range of applicability of 
second order viscous hydrodynamics, corrections from bulk viscosity are found to be 
small compared to those from shear viscosity if the bulk relaxation time peaks around 
Tc. This means that if bulk and shear viscous correction do not become so large as to 
render the hydrodynamic expansion invalid for the relevant part of a heavy-ion collision, 
the extraction of the shear viscosity should be possible to reasonable accuracy when 
neglecting the effect of bulk viscosity. 

Apart from viscous corrections to the thermal equilibrium distribution functions 
becoming large, bulk and shear viscous corrections can lead to a negative longitudinal 
pressure, and hence the breakup of the system into droplets [311 ESI EH]. Second order 
viscous hydrodynamics can hint at when such cavitation happens but is not suited to 
describe the process as it happens outside the range of its applicability. 

5. Progress beyond second order viscous hydrodynamics 

The problem with the viscous correction to the stress energy tensor becoming larger 
than the equilibrium part, leading to negative pressure, arises especially in the early 
stage of the evolution, when the local momentum distribution is not yet equilibrated 
but highly anisotropic due to rapid longitudinal expansion. Recent progress has been 
made in describing this early time evolution and late time hydrodynamics within the 
same framework by performing the hydrodynamic expansion around an anisotropic 
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distribution [371 [Ml ESI SO]- This results in new equations of motion, including one 
for the degree of anisotropy of the distribution function. The procedure can reproduce 
both the limits of free streaming and ideal hydrodynamics and results in second order 
viscous hydrodynamics when expanding around a small anisotropy parameter, which so 
far has been shown in the one dimensional case |37j . 

Another logical step is to expand to third order in gradients as done in |¥Tj . 
Numerical differences to second order Israel Stewart theory are completely negligible 
for Tj/s = 0.05 and become significant for rj/s > 0.2. It has also been pointed out that 
the hydrodynamic equations depend on the details of their derivation. While Israel and 
Stewart used the second moment of the Boltzmann equation to derive hydrodynamic 
equations for the dissipative currents, in [32] the definition of the latter was used directly. 
This leads to equations of motion of the same form but with different coefficients. 
Microscopic transport calculations [^ Hlj show very good agreement with the new 
equations of motion up to rj/s ~ 3, while the Israel-Stewart equations show large 
differences for rj/s > 0.2. This demonstrates that the details of the derivation are 
relevant in particular when hydrodynamics is being matched to kinetic theory at late 
times. 

6. Temperature dependent ri/s 

What is extracted from experimental data by comparison with viscous hydrodynamic 
calculations using constant rj/s is at best an effective, or average, {rj/s). In reality r}/s 
should depend on the local temperature of the medium, dropping from large values at 
high temperatures to a minimum at T^, and rising with decreasing temperature in the 
hadronic phase [45ll46ll47l l48] . 

It has been shown that when including such modeled temperature dependence 
in calculations at RHIC energies, the details of ri/s{T) in the quark-gluon plasma 
phase have little influence on the final elliptic flow result, while hadronic ri/s{T) 
modifies V2 strongly [H]. At the highest LHC energies the conclusion is the opposite: 
weak dependence on the hadronic r]/s but strong dependence on r]/s{T) in the QGP 
phase. Interestingly, at RHIC energies, a significant dependence has been found on the 
minimum value of rj/ s{T) around [l9]. This might indicate that we are determining 
such a minimal value, rather than (rj/s), when using a constant rj/s. 

In addition, it was pointed out in [50] that there is a strong dependence on the 
initial value of n^'^ when starting with a large rj/s. So, particularly when including a 
temperature dependent rj/s, it is essential to gain a better understanding of the pre- 
equilibrium stage in heavy-ion collisions to determine the initial conditions for viscous 
hydrodynamics. 
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7. Viscous corrections to particle distribution functions 

When translating the dissipative T'^'^ to particles in the Cooper- Frye formalism [5T] . 
corrections to the distribution function 6f have to be taken into account: 

for an component system with dn the degeneracy of species n. The ansatz 

5/n=^/o(l±/or/x(p)^, (8) 

where is a unit vector in the a direction, leaves some freedom and the usual procedure 
is to assume that all coefficients C„ are equal, even though they should depend on the 
individual particle species' interaction rate (see [52]), and use x{p) = which is derived 
within a relaxation time plus Boltzmann approximation: 

5fn = /on(l ± fon)py^7r^0^^^^:py^ V n . (9) 

It has been shown in [53] that x{p) where a can take on values from 1 to 

2, which is the case for example for a system with radiative and elastic energy loss 
that has x{p) pi-^s. xhe exact form of the correction has a large effect on the pt 
differential elliptic flow for Pt ^ 1 GeV. It should be noted, however, that the analysis 
of experimental data in [51] indicates x{p) p"^- A more general problem is that 
corrections 6f can become large compared to /q. For hadrons this problem can be 
reduced when using a hadronic afterburner and switching at intermediate temperatures 
of ~ 160 MeV, but for photons that are produced throughout the whole evolution this 
becomes a serious concern [231. 



8. Hadronic afterburner 



The use of hybrid models coupling early hydrodynamic evolution to a microscopic 
hadronic cascade in the later stage has a long history [551 EHl EH [591 l60l [61] . Recently 
2+1 dimensional viscous hydrodynamics has been coupled to UrQMD [621I631 IM I I65 1 I66] . 
When using a chemically frozen equation of state, transverse momentum spectra of 
produced particles show very little dependence on the temperature Ts„ at which 
the switching between hydro and UrQMD is performed. However, V2 does show a 
dependence unless an increasing t]/s{T) is introduced in the hadronic stage of the hydro 
simulation to emulate the effect of UrQMD's larger dissipation [62]. This rj/s rises 
to ~ 0.4, for which the applicability of second order viscous hydrodynamics becomes 
questionable. Also, it turns out that the dynamics in UrQMD cannot be described by 
viscous hydrodynamics with r]/s{T) due to different relaxation times [62], underlining 
the need for the more realistic microscopic rescattering. 
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9. Initial conditions and event-by-event hydrodynamics 

Fluctuating initial conditions for hydrodynamic simulations of heavy-ion collisions 
have been argued to be very important for the exact determination of collective flow 
observables and to describe features of multi-particle correlation measurements in heavy- 
ion collisions [nilEHlEHllTolinillSlITallHIISlITSlI^ 

Real event-by-event hydrodynamic simulations have been performed and show 
modifications to spectra and flow from "single-shot" hydrodynamics with averaged 
initial conditions [8T| [T6| \T7\ [8^ [85] . An important advantage of event-by-event 
hydrodynamic calculations is the possibility to study higher flow harmonics such as 
^3, which are entirely due to fluctuations in the initial conditions. Different f„ depend 
differently on 77/5 and the details of the initial condition, which is determined by the 
dynamics and fluctuations of partons in the incoming nuclear wave functions. This 
observation can be used to determine these long sought after details of the initial state 
and medium properties in heavy-ion collisions by performing a systematic analysis of all 
harmonics up to e.g. n = 6 as a function of rj/s and the initial state properties and 
compare to experimental data. First predictions of V3 [16] agree extremely well with 
experimental data from RHIC [HE]. Furthermore, it has now been shown that at low pt 
(and |A?7| > 1 (ALICE), |A?7| > 2 (ATLAS)), the main features of dihadron correlations 
in the angular difference A0 between the hadron momenta can be described by flow, i.e. 
the sum of Vi to Vq only [HTJ [88]. The double-peak structure on the away-side is hence 
described mostly by (triangular) flow as predicted in [771 179] . 

10. Summary of LHC predictions, conclusions and outlook 

Calculations using ideal and viscous hydrodynamics at LHC energies hint at little or no 
increase of (rj/s) [HH [90], [6OI [611 [TTl [651 [SI]- However, to be consistent a temperature 
dependent rj/s should be employed to extend RHIC calculations to LHC (see e.g. [65]). 
Also, essentially all models underestimate f2(pr) for px ^ 0.8 GeV, which might be 
explained by contributions from non-thermalized particles [92]. The issue is however 
not settled yet. Nevertheless, the overall early success of hydrodynamics in describing 
the main features of flow of charged hadrons at the LHC indicates that the QGP at LHC 
is also a nearly perfect fluid. Higher harmonics are also being computed for the LHC 
and comparison of experimental data to flrst predictions from event-by-event simulations 
[T7] show that hydrodynamic results are in the right ball park. As mentioned earlier, 
a systematic analysis will be needed to potentially determine ri/s{T) and the details of 
the initial conditions at LHC with the use of all 

Better understanding of viscous corrections to the distribution functions and 
dependence on the model for coupling to a hadronic afterburner, detailed studies of 
3-1-1 dimensional hydrodynamics with viscosity, and above all a better understanding of 
the pre-equilibrium stage and its transition to hydrodynamics are the next important 
steps in developing a reliable viscous hydrodynamic description of the bulk dynamics in 
heavy-ion collisions. 
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